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1 Introduction 

The minimal supersymmetric standard model (MSSM), like the standard model itself, pre- 
dicts a zero mass for neutrinos, and this is not compatible with the recent neutrino ob- 
servations. One of the most promising methods to attribute a tiny but nonzero mass to 
neutrinos is the seesaw mechanism [1], which requires three extremely heavy right-handed 
neutrinos. In the MSSM with right-handed neutrino supermultiplets, the leptonic part of 
the superpotential is 

W = Y; j e a(3 H?l c Ri L? + Y^e ap H^v Ri L^ + l -M l0 v Ri v Rj , (1) 
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where Lj is the supermultiplet corresponding to the doublet (vLj, lij)- Without loss of 
generality, we can rephase and rotate the fields to make the matrices Y{ 3 and Mij real and 
diagonal: Y{ 3 = diag(y e , Y^, Y T ) and M* J = diag(Mi, M2, M3). In this basis, Y v can have 
off-diagonal and complex elements. Soft supersymmetry breaking terms of the Lagrangian 
in the context of this model can include 

= (mlY^ui + [\B V M^^ R + H.c.]. (2) 

The second term in Eq. (2), the "neutrino 5-term," is a lepton- number- violating term [2] 
which can cause profound effects including sneutrino-antisneutrino oscillation [3, 4]. The 
parameter B v is allowed to be much larger than the electroweak scale because it is associated 
only with vr, which is an electroweak singlet. It has been shown that a nonzero neutrino 
S-term can create neutrino mass through one-loop diagrams [3]. The upper bound on the 
neutrino mass can then be translated into a bound on B u , 

B u < 10 3 m . (3) 

If B v is large, some new effects are expected both in the e + e~ accelerator experiments [3] 
and in cosmology [4] . In particular, values of B u close to the saturating bound (3) can induce 
observable slepton-antislepton oscillation. In this paper, we show that large values of B u can 
also affect other observables. 

It is well-known that nonzero flavor-number-violating slepton mass terms in the soft 
Lagrangian (m^gL^L^, a 7^ (3) can give rise to rare decays such as (/i — > je), (r — > je), 
and (r — > 7/1) . One way to avoid flavor changing neutral current (FCNC) effects is to choose 
the off-diagonal mass terms to be small. In fact, theories such as minimal supergravity 
(mSUGRA) suggest that at the grand unified theory (GUT) scale, the soft supersymmetry 
breaking terms are flavor blind; that is, at the GUT scale 

-C^ = ml(LL x L La + Fj: i to + i>l x i> I to + 

+ {bH x H 2 + H.c.) + a Q {Yi 3 e a3 Hfl^L% + Y?e a3 H«u m L%) 

+ Q^M^ + H.c), (4) 

with universal mi/2, and do- 

The off-diagonal elements of the neutrino Yukawa coupling radiatively produce nonvan- 
ishing off-diagonal mass terms for the left-handed slepton doublet: 

Y ka {Y kPy j- / r AcJJT ^ \ rA Gl/T i 2 l 
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This effect was first noticed and studied in [5] and then worked out in a series of papers 
(e.g., see [6]). However, the contribution to m^p from the neutrino 5-term was ignored. In 
section 2, we study this effect and show that if B v is large, its contribution will dominate 
over the effects in Eq. (5). 

In the MSSM with flavor blind soft supersymmetry breaking terms, in addition to the 
phases in the Yukawa couplings, there are two other independent CP- violating phases, usually 
chosen to be the phases of the ao and fi parameters. These phases can create electric 
dipole moments (EDMs) for charged leptons and for the neutron [7]. In the presence of the 
neutrino i?-term, there is one more phase which can also give a contribution to the EDM of 
charged leptons. In section 3, we show that, even if at the GUT scale no A-term is present 
(ao = 0), through one-loop corrections, the neutrino 5-term creates A-terms for leptons at 
the electroweak scale. This effect could be the dominant term in lepton EDMs. 

In section 4, we explore the upper bounds on the imaginary and real parts of B v resulting 
from the upper bounds on the branching ratios of the rare decays [BR(Z a — > Z 7 + 7)] and the 
EDMs of the charged leptons. The main limitation will be the uncertainty in the pattern of 
neutrino Yukawa couplings Y v . 

2 Effects of the neutrino B-term on slepton mixing 

It has been shown that the off-diagonal slepton masses (m^L^L^, a 7^ (3) at the one-loop 
level can give rise to lepton-number- violating rare decays such as (fx — > ej), (r — > /ry), and 
(t — > ej) [6, 8]. In the mass insertion approximation, a simplified formula can be derived 
[9]: 

Br(l a - l p + 7) ~ ^^^tan 2 /3. (6) 

m susy 

The upper bounds on the branching ratios of the rare decays [10] can be interpreted as 
bounds on the off-diagonal elements of |77i^|: 

I 2 I <r- ^ X ( msus D \2 2 I 2 1 0-4 . Tn susy 2 2 (rj\ 

|m ^' tan/5 l 200 GeV J ^ ta^ l 200 GeV J msusy [ ' 



and 



I 2 I > 1 / m susy \2 2 /o\ 

|m ™ 1 < W ( 200-GeV } (8) 



The next generation of experiments [11] is expected to improve Eq. (7) to 



I 2 I ^ 6 X 10 ITlsusy \2 2 i 2 | 0-07 , fn susy 2 2 /q\ 

tan/3 l 200 GeV J ^ ta^ l 200 GeV J 1 J 



3 



The off-diagonal mass terms for left-handed sleptons receive a contribution from the 
neutrino B-term through the two diagrams shown in Fig. 1. The two lepton number violating 
vertices on the neutrino line are the neutrino B-term and the standard vr mass term. The 
neutrino A-term is also needed. The amplitude corresponding to diagram (a) is equal to 

d A k 



iM k i 



(4tt) 

Similarly, diagram (b) gives 



E k — 2 alY^{Y^yB v . (10) 



The mass correction is given by the sum of the two amplitudes: 



Y u ka (Y^yRe[a B;] 
™(2)a/3 = "2 L > ( 12 ) 



which has to be added to m 2 ^^ presented in Eq. (5). For Re[a -B*] ~ 10ml, rn 2 2 } a p exceeds 
m \i)ap- Note that the contribution we have found does not depend on the heavy right- 
handed masses at all. This can be traced back to the form of the neutrino .B-term assumed 
in Eq. (4). Had we defined this term as B 2 vrVr, the result would have been proportional 

to Re[a (B 2 u )*]/M k . 

Up to factors of \og(A G uT / M k ) , m 2 i) a p and m 2 2 ^ a/3 [see Eqs. (5,12)] have the same flavor 
structure. The structure can be different only if the masses of right-handed neutrinos are 
hierarchical (M 1 <C M 2 <C M 3 ). Although the one- loop mass matrix presented in Eq. (5) is 
enhanced by a factor of 6 \og(A GUT /M k ) ~ 10, the neutrino S-term contributions given in 
Eq. (12) dominate if B v ~ 10 3 m as allowed by Eq. (3). 

The dependence of ra 2 2 ^ al3 on B u involves the combinations J2k Y^iY^)*. To derive 
conclusive bounds on B u , first we have to find some lower bounds on the J2k Y^ a {Y^)* 
combinations; however, this information is not available at the moment. If Y u are so large 
that m 2 i) al 3 [see Eq. (5)] saturate the bounds (7,8), B v has to be smaller than 10m [to be 
compared with Eq. (3)]. In section 4, we will discuss this further. 

In the discussion above, we have assumed a ~ m ~ m susy but it is possible that a 
is much smaller than m . In this case, m 2 2 ^ a/3 given in Eq. (12) will not be the dominant 
effect. At the one loop level, there is no contribution to m 2 al3 proportional to \B U \ 2 : it can 
be shown that the two one-loop diagrams that are proportional to \B V \ 2 (depicted in Fig. 2) 
cancel each other at zero external momentum. However, at the two-loop level, there is a 
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contribution proportional to \B U \ 2 which can dominate over fn^ a/3 [Eq. (5)] provided that 

\b v \ 2 y u y;i^) 2 > m 2 . 

3 Effects of the neutrino B-term on A-terms 

In this section, we show that the neutrino .B-term creates an A-term for leptons through 
one-loop diagrams. We then discuss how this will affect the EDMs. 

As is discussed in the literature [7], the phases of // and a$ can create electric dipole 
moments for charged leptons as well as for the neutron. The current bounds on lepton 
EDMs are 

d e < 1.5 x 1(T 27 e cm [12] < 7 x 1(T 19 e cm [10] (13) 

and 

d T <3x 1(T 16 e cm [10]. (14) 
Proposed future experiments are expected to set stronger bounds: 

d e < 10~ 32 e cm [13] d^ < 10~ 24 (5 x 10~ 26 ) e cm [14] ([15]). (15) 

The bounds on the electric dipole moments of the charged leptons yield strong bounds on 
the imaginary parts of /i and ao [16]. 

The phase of the neutrino B-term can provide yet another source of CP-violation, t 
When we fixed the mass matrix M to be real, the phases of vr were fixed; therefore, the 
phase of B u in this convention cannot be removed. We expect the imaginary part of B v to 
give contribution to EDMs. 

The parameter B v contributes to the A^-term through the diagram shown in Fig. 3. 
Adding this correction to the tree level A\ [see Eq. (4)], we find 

= -iaoY',% - ^yi« ( y«)-i;«B„. (16) 

Similarly, the neutrino .B-term contributes to the the A u -teim through the diagram shown 
in Fig. 4: 



q 



^Within the extended MSSM, the Yukawa couplings (Y v ) are another source of CP-violation. This effect 
has been discussed in [17]. 
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= -iaoY?-j^Y?(Y?yY*B v . (17) 

According to [16], for m ~ 200 GeV, the present bound (d e < 10 -27 e cm) implies 

Im(in(ffi>/Km )^0.1. (18) 

Since the dependence of d e on lm(Af e ) is linear, if in the future the bound d e < 10~ 32 e cm 
is obtained, the above bound will be improved to 

Im(An(^i)/Km )~10- 6 . (19) 

Assuming that B u is the only source of CP- violation, the bound in Eq. (18) can be translated 
into 

Im^] J2 Yr{Y v k yY u ke < 15m m e /(iJ 1 ), (20) 

k 

which can be improved by five orders of magnitude in the future. 

The present experimental data do not lead to any conclusive bounds on the values of 
the neutrino Yukawa couplings. In section 4, we will discuss how future observations and 
developments can improve our knowledge on Y v . In principle, Y v can be as large as order 1. 
(In fact, Y u can be even larger than 1; however in this case we cannot treat it perturbatively.) 
For Y ke (Y ke )* ~ 1, the present bound (18) gives Im[5„] < 10m [we have used Eq. (16) and 
m e = (i?i)yj ee ]. Future EDM experiments can make the bound dramatically stronger. 

Discovery of a lepton EDM could provide invaluable information on B u . The neutrino 
5-term gives a contribution to Af, however, it has no impact on the A-term of quarks. As 
a result, lm(B u ) will not affect the EDM of the neutron. On the other hand, Im(a ) and 
Im(/x) give contributions to both the EDM of charged leptons and neutron. It is possible 
that the contributions of Im(a ) and Im(/i) cancel each other. However, it has been shown 
[18] that cancellation in the electron EDM occurs in the same regions as cancellation in the 
neutron EDM. Therefore, if d e turns out to be nonzero while d n <C d e , the effect cannot be 
attributed to the contribution of Im(a ) or Im(/x). Such a situation can be explained with a 
nonzero complex B v . 

There is another point that is noteworthy: EDMs are proportional to Y J k\Y ka \ 2 - l and 
BR(/ a -> Ipj) are given by \E k Y ka (Y^)*\ 2 , which are both independent of M h To the 
author's best knowledge, there is no other observable that depends on these combinations. 
If \B V \ is large, by studying these observables we can extract additional information on 
Yukawa couplings which will improve our current understanding of the seesaw mechanism 
and leptogenesis. 
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If the neutrino B-term gives the dominant contribution to the electric dipole moments, 
we expect d T /(m T £ fc \Y* T \ 2 ) = d M /(m M £ fe |Y^| 2 ) = d e /(m e £ fe |^ fee | 2 ); therefore, if d e is 
close to its present upper bound, d e ~ 10 -27 e cm, we expect d^ ~ 10~ 25 e cm which can be 
tested in proposed experiments [15]. 

4 Bounds on B u 

In sections 2 and 3, we have shown that large values of B v can lead to flavor-violating rare 
decays and EDMs of charged leptons. However, the dependence of these observables on B v is 
through the unknown combination of Yukawa couplings Y„(Y*P)*. To derive upper bounds 
on B u , we have to find other observables that provide lower bounds on these combinations. 
In this section, we combine various pieces of information on the Yukawa couplings (some of 
them yet to be obtained) to derive a lower bound on the factors Y* a (Y*P)*. We will then 
use the current upper bounds on the branching ratios of the rare decays and the values of 
EDMs to extract upper bounds on B u . 

Neutrino masses depend on the Yukawa couplings through 

™$ = En te ^n*W. (21) 
k ivik 

All the parameters involved have to be evaluated at the electroweak scale. The effect of 
running from the GUT scale to the electroweak scale can change the details of the neutrino 
masses and mixing [19], however, the order of magnitude of the masses will not be affected. 
Since, here, we want to estimate only the order of magnitude of the Yukawa couplings, we 
will neglect the running effects. 

Currently we have only bounds on the neutrino masses [20]: 

\j&rn 2 atm < m v < 1 eV, 

where Am 2 tm = 2.5 x 10~ 3 eV 2 [21]. Future terrestrial and cosmological experiments will 
improve these bounds. Our knowledge of the masses of the right-handed neutrinos (Mfc) is 
even less complete than the information on m^p- If leptogenesis is the mechanism behind 
the baryon asymmetry of the universe [22] , it will be possible to derive a lower bound on M k 
[23]. For a given value of m^p, there is at least one k such that 

\Y^(Y^Y\ > ± ( m * ) 2 X 10 ' 6 ( Mk ) (22) 
1 » [ " } 1 > 3 ^0.1 eV) sin 2 (3 U x 10» GeV J " 1 ] 



7 



The parameter m^) can be extracted directly from neutrinoless double decay observa- 
tions. If m^) is relatively large (m^) ~ 0.1 eV), its effect should be observable in future 
experiments [24]. Using Eq. (16), d e < 1CT 27 e cm [see Eq. (18)] yields 

M^) <3 xl0'(^) { 6 X 108 GeV ) . 
m \ mil? A M " I 

In future, if the bound d e < 1CT 32 is obtained, this bound will be improved by five orders of 
magnitude [see Eq. (19)] which means it can be more restrictive than the bound in Eq. (3). 

Extracting bounds on Y^ a (Y^)* , a ^ f3, will be more challenging because, unlike m'fj, 
m^p, a ^ [3, cannot be directly measured. Even if forthcoming experiments find that the 
overall neutrino mass is of order of 0.1 eV (quasidegenerate mass scheme), it will be difficult 
to derive definite lower bounds on [m^j, (a ^ f3)\. In the case of a quasidegenerate mass 
scheme (the scheme for which the absolute values of the mass eigenvalues are much larger 
than \J Am 2 atm ) with zero Majorana phases we expect m^,m^),m^) <C m^J , mfy , m£} . 
Only in the framework of the quasidegenerate neutrino mass scheme with at least one nonzero 
Majorana phase, it is possible to have large off-diagonal neutrino masses, mW,m^),mfg) 3> 
\J Aml tm (a phase equal to n also works). On the other hand, determining the values of 
Majorana phases is very challenging, if possible at all [25]. Nevertheless, let us suppose that 
in the future some hypothetical experiment will be able to determine m^}. Then, assuming 
that the factors Y^ e (Y^)* do not cancel each other, the present bound on \m\ | given in 
Eq. (7) implies 

or D *w 2 10 5 / m \ 2 /(U eV\ /6 x 10 8 GeV\ /OQ , 
Re[a B;}/ml < — — — — —j-- . (23) 



tan/? V200 GeV/ \ mH$ I \ M i 



k 



The future possible bounds [inferred from Eq. (9)] can be more restrictive than Eq. (3): 

t~> r d*i / 2 3xl0 3 / m \ 2 f 0.1 eV\ /6 x 10 8 GeV\ , nA . 

however, as we pointed out earlier, at the moment, measuring m'fj seems to be impossible. 

5 Concluding remarks 

We have studied the effects of the neutrino 5-term on the slepton mixing and EDMs of 
charged leptons in the framework of seesaw model embedded in the MSSM with universal 
soft supersymmetry breaking terms. 

If B u > 10m ~ 10a but a > B U Y U (Y U )* /(An) 2 , the dominant flavor-violating slepton 
masses are given by Eq. (12) rather than Eq. (5). For values of B v satisfying B V Y V {Y V )* / (An) 2 > 
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do and \B V \ 2 Y V {Y V )* / (4-7r) 2 > itlq, the two-loop contribution proportional to \B U \ 2 can be dom- 
inant. The bounds on the Yukawa couplings and neutrino S-term which have been discussed 
in the literature allow quite large contributions to the slepton masses, violating the upper 
bounds from rare flavor-violating decays. However, since there is no direct lower bound 
on the combinations of neutrino Yukawa couplings appearing in the formulations, it is not 
possible to derive any upper bound on the Re[a -B*]. 

The parameter B v can be considered as another source for CP-violation and therefore 
EDMs. In fact, we have shown that the neutrino 5-term directly creates A-terms both 
for neutrinos and charged leptons-but not quarks-even if = at the GUT scale. The 
imaginary part of A\ gives a contribution to the EDMs of charged leptons. If the B v effect is 
dominant, we expect d T /(m T J2k l^ fcT | 2 ) — <W("fyZ)fc |^ fc/1 | 2 ) = d e /(m e J2k l^ fce | 2 ); therefore, 
if d e is close to its present upper bound we expect that the proposed experiments [15] will 
be able to measure the value of d^. The discovery of nonzero d e and d^ while d n <C d e 
can be explained with large \m(B u ). In this case, d e oc J2k Y^ e (Y^ e )*lm(B l/ ). If lm(B u ) is 
determined through some other observation, information on d e and lm(B v ) combined with 
m^ee = Y J k(H2) 2 (Yl; e ) 2 / M k (extracted from neutrinoless double beta decay searches) can 
provide us with information on the values of M k , shedding light on the origins of neutrino 
masses and on leptogenesis. 

Acknowledgment 

I would like to thank Y. Grossman, M. Peskin, and A. Yu. Smirnov for useful discussions 
and encouragement. I am also grateful to M. M. Sheikh- Jabbari for careful reading of the 
manuscript. 

References 

[1] T. Yanagida, in Proceeding of Workshop on Unified Theory and Baryon Number of the 
Universe, Tsukuba, Japan, 1979, edited by O. Sawada and A. Sugamoto (KEK, 1979) 
p. 95; M. Gell-Mann, P. Ramond and R. Slansky, in Supergravity, edited by P. van 
Niewwnhuizen and D. Freedman (North Holland, Amsterdam, 1979); R. N. Mohapatra 
and G. Senjanovic, Phys. Rev. Lett. 44 (1980) 912. 

[2] M. Hirsch, H. V. Klapdor-Kleingrothaus, and S. G. Kovalenko, Phys. Lett. B 398 
(1997) 311; Phys. Lett. B 403 (1997) 291; 



9 



[3] Y. Grossman and H. E. Haber, Phys. Rev. Lett. 78 (1997) 3438. 

[4] L. J. Hall, T. Moroi, and H. Murayama, Phys. Lett. B 424 (1998) 305; Y. Grossman, 
T. Kashti, Y. Nir, and E. Roulet, Phys. Rev. Lett. 91 (2003) 251801; G. D'Ambrosio, 
G. F. Giudice, and M. Raidal, Phys. Lett. B 575 (2003) 75. 

[5] F. Borzumati and A. Masiero, Phys. Rev. Lett. 57 (1986) 961. 

[6] J. Hisano, T. Moroi, K. Tobe and M. Yamaguchi, Phys. Rev. D 53 (1996) 2442. 

[7] T. Moroi, Phys. Rev. D 53 (1996) 6565; Erratum-ibid. D 56 (1997) 4424; S. Pokorski, 
J. Rosiek and C. A. Savoy, Nucl. Phys. B 570 (2000) 81; J. L. Feng, K. T. Matchev and 
Y. Shadmi, Nucl. Phys. B 613 (2001) 366; M. Garena, et al, Phys. Lett.B 390 (1997) 
234; E. Gabrielli and U. Sadrid, Phys. Rev. Lett. 79 (1997) 4752; J. L. Feng and K. T. 
Matchev, Phys. Rev. Lett. 86 (2001) 3480; L. Everett et al, Phys. Rev. Lett. 86 (2001) 
3484; T. Ibrahim et al, Phys. Rev.D 64 (2001) 016010; J. Ellis et al, Phys. Lett.B 507 

(2001) 224; Z. Chacko and G. D. Kribs Phys. Rev.D 64 (2001) 75015; D. G. Gerdeno et 
al, Phys. Rev.D 64 (2001) 093012; U. Chattopadhyay and P. Nath, Phys. Rev. D 66 

(2002) 093001; S. P. Martin and J. D. Wells, Phys. Rev. D 67 (2003) 015002; T. Ibrahim 
and P. Nath, Phys. Rev.D 58 (1998) 111301; T. Falk and K. Olive, Phys. Lett.B 439 
(1998) 71; M. Brhlik, G. Good and G. L. Kane, Phys. Rev.D 59 (1999) 115004; U. 
Chattopadhyay, T. Ibrahim and P. Roy, Phys. Rev. D 64 (2001) 013004; V. Barger et 
al, Phys. Rev. D 64 (2001) 056007; S. Abel, S. Khalil and O. Lebedev, Nucl. Phys. 
B 606 (2001) 151; T. Ibrahim and P. Nath, Phys. Rev. 64 (2001) 093002. 

[8] J. Hisano and D. Nomura, Phys. Rev. D 59 (1999) 116005. 

[9] J. A. Casas and A. Ibarra, hep-ph/0109161; 

[10] K. Hagiwara et al, Phys. Rev. D 66 (2002) 010001; K. Abe et al. [Belle Collaboration], 
arXiv:hep-ex/0310029. 

[11] T. Oshima, Prepared for 3rd Workshop on Neutrino Oscillations and Their Origin 
(NOON 2001), Kashiwa, Japan, 5-8 Dec 2001; L. M. Barkov et al, proposal for an 
experiment at PSI, http://meg.web.psi.ch. I. Hinchliffe and F. E. Paige, Phys Rev. 
D 63 (2001) 115006; D. F. Carvalho et al, hep-ph/0206148; J. Kalinowski, Acta Phys. 
Polon. B 33 (2002) 2613. 

[12] E. D. Commins et al, Phys. Rev. A 50 (1994) 2960. 



10 



[13] S. K. Lamoreaux, nucl-ex/0109014. 

[14] R. Garey et al, letter of intent to BNL (2000); Y. K. Semertzidis et al, hep-ph/0012087. 

[15] J. Aysto et al, hep-ph/0109217. 

[16] I. Masina and C. A. Savoy Nucl. Phys. B 661 (2003) 365 [arXiv:hep-ph/0211283]. 

[17] J. R. Ellis, J. Hisano, S. Lola and M. Raidal, Nucl. Phys. B 621, 208 (2002) 
[arXiv:hep-ph/0109125]; J. R. Ellis, J. Hisano, M. Raidal and Y. Shimizu, Phys. Lett. 
B 528, 86 (2002) [arXiv:hep-ph/01 11324]; Phys. Rev. D 66, 115013 (2002) [arXiv:hep- 
ph/0206110]; J. R. Ellis and M. Raidal, Nucl. Phys. B 643, 229 (2002) [arXiv:hep- 
ph/0206174]; I. Masina, Nucl. Phys. B 671, 432 (2003) [arXiv:hep-ph/0304299]; Y. 
Farzan and M. E. Peskin, work in preparation. 

[18] T. Falk, K. A. Olive, M. Pospelov, and R. Roiban, Nucl. Phys. B 560, 3 (1999) 
[arXiv:hep-ph/9904393]. 

[19] M. Frigerio and A. Y. Smirnov, JHEP 0302, 004 (2003) [arXiv:hep-ph/0212263] and 
references therein. 

[20] S. Hannestad, JCAP 0305:004,2003 [astro-ph/0303076]; J. Bonn et al, Nucl. Phys. 
Proc. Suppl. 91 (2001) 273; V. M. Lobashev et al, Proc. of the Int. Con. Neutrino 
2000, Sudbury, Canada; Nucl. Proc. Suppl. 77 (1999) 327; Nucl Phys. Proc. Suppl. 91 
(2001) 280; V. M. Lobashev, Prog. Part. Phys. 48 (2002)123. 

[21] Super-Kamiokande Collaboration, T. Nakaya, eConf C020620 (2002) SAAT01, [hep- 
ex/0209036]. 

[22] M. Fukugita and T. Yanagida, Phys. Lett. B 174 (1986) 45. 

[23] E. Kh. Akhmedov, M. Frigerio and A. Yu. Smirnov, JHEP 0309:021,2003; W. Buch- 
muller, P. Di Bari and M. Plumacher, Nucl. Phys. B 643 (2002) 367; S. Davidson and 
A. Ibarra, Phys. Lett. B 535 (2002) 25. 

[24] S. Pirro et al, Nucl Instrum. Meth. A444 (2000) 71; E. Fiorini et al, Phys. Rep. 307 
(1998) 309; A. Bettini, Nucl. Phys. Proc. Suppl. 100 (2001) 332; H. Ejiri, J. Engel, H. 
Hazama, P. Krastev, N. Kumodi and R. G. H. Robertson, Phys. Rev. Lett. 85 (2000) 
2917; GENIUS Collaboration, H. V. Klapdor-Kleingrothaus et al, hep-ph/9910205. 



11 



[25] V. Barger, S.L. Glashow, P. Langacker, and D. Marfatia, Phys. Lett. B 540 (2002) 
247; H. Nunokawa, W.J.C. Teves and R. Z. Funchal Phys. Rev. D 66 (2002) 093010; 
F. Feruglio, A. Strumia and F. Vissani, Nucl. Phys. B 637 (2002) 345. 



12 





13 



H 2 




Figure 3: Diagram contributing to A t . and F[ represent the auxiliary fields associated 
with v\ and Lj, respectively. 




Figure 4: Diagram contributing to A u . F£ and F[ represent the auxiliary fields associated 
with v q R and Lj, respectively. 
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